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We search for evidence of a light scalar boson in the radiative decays of the Y{2S) and T{3S) 
resonances: T{2S,3S) — > 7^", A" fi'^ fi~ . Such a particle appears in extensions of the Standard 
Model, where a light CP-odd Higgs boson naturally couples strongly to 6-quarks. Wo find no 
evidence for such processes in the mass range 0.212 < m^o < 9.3 GcV in the samples of 99 x 10*' 
T{2S) and 122 x 10*^ r{3S) decays coUected by the BABAR detector at the PEP-II B-factory and set 
stringent upper limits on the effective coupling of the b quark to the A°. We also limit the dimuon 
branching fraction of the r]b meson: B{r]b — > It^ lJ-~) < 0.9% at 90% confidence level. 

PACS numbers: 13.20.Gd, 14.40.Gx, 14.80.Cp, 14.80.Mz, 12.60.Pr, 12.15.Ji 



The concept of mass is one of the most intuitive ideas ence. Yet the fundamental nature of mass remains one 
in physics since it is present in everyday human experi- of the great mysteries of science. The Higgs mechanism 
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is a theoretically appealing way to account for the dif- 
ferent masses of elementary particles It implies the 
existence of at least one new scalar particle, the Higgs bo- 
son, which is the only Standard Model (SM) particle 
yet to be observed. The SM Higgs boson mass is con- 
strained to be of 0(100 - 200 GeV) by direct searches 0] 
and by precision electroweak measurements 

A number of theoretical models extend the Higgs sector 
to include additional Higgs fields, some of them naturally 
light 5]. Similar light scalar states, e.g. axions, appear 
in models motivated by astrophysical observations and 
are typically assumed to have Higgs-like couplings Q. 
Direct searches typically constrain the mass of such a 
light particle, to be below 2mh Q, making it ac- 
cessible to radiative decays of T resonances Q. Model 
predictions for the branching fraction (BF) of T ^ 7^4° 
decays range from 10~^ @, 0] to as high as 10""* @. 
Empirical motivation for a low-mass Higgs search comes 
from the HyperCP experiment [lol | , which observed three 
anomalous events in the 17+ pjX^ ii~ final state. These 
events have been interpreted as production of a scalar bo- 
son with the mass of 214.3 MeV decaying into a pair of 
The large datasets available at BaBAR allow 
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muons 

us to place stringent constraints on such models. 

If a light scalar exists, the pattern of its decays 
would depend on its mass. Assuming no invisible (neu- 
tralino) decays [l3|, for low masses m^o < 2tot- the BF 
B^fj, = B{A^ should be sizable. Significantly 

above the r threshold, A" t^t~ would dominate, and 
hadronic decays might also be significant. 

This Letter describes a search for a resonance in the 
dimuon invariant mass distribution for the fully recon- 
structed final state T{2S,3S) 7^°, ^ ^l+fl-. We 
assume that the decay width of the resonance is negligi- 
bly small compared with the experimental resolution, as 
expect ed Ig , Il3j for m^o sufficiently far from the mass of 
the [141. We further assume that the resonance is a 
scalar (or pseudo-scalar) particle. While the significance 
of any observation would not depend on this assumption, 
the signal efficiency and, therefore, the BFs are computed 
for a spin-O particle. In addition, following the recent dis- 
covery of the r]b meson [31 , we look for the leptonic decay 
of the rjh through T{2S, 3S) ^r]i,, -qb We use 

T{r]b) = 10 ± 5 MeV, the range expected in most theoret- 
ical models and consistent with the BaBAR results [3l • 

We search for two-body transitions T{2S, 35') 7^°, 
followed by decay ^-^ in samples of (98.6 ± 

0.9) X 10^ T{2S) and (121.8 ± 1.2) x 10^ T{3S) de- 
cays collected with the BABMi detector at the PEP-II 
asymmetric-energy e+e" collider at the SLAG National 
Accelerator Laboratory. We use a sample of 79fb~^ ac- 
cumulated on the T(45') resonance (T(4S') sample) for 
studies of the continuum backgrounds. Since the T(4S') 
is three orders of magnitude broader than the T{2S) and 
YiiS), the BF B{T{AS) 7^°) is expected to be negli- 
gible. For characterization of the background events and 



selection optimization, we also use a sample of 1.4fb ^ 
(2.4fb"^) collected 30 MeV below the T{2S) {Y{iS)) res- 
onance (off-resonance samples). The BABAR detector is 
described in detail elsewhere [3, [3] ■ 

We select events with exactly two oppositely-charged 
tracks and a single energetic photon with a center-of- 
mass (CM) energy E* > 0.2 GeV, while allowing addi- 
tional photons with GM energies below 0.2 GeV to be 
present in the event. We assign a muon mass hypoth- 
esis to the two tracks (henceforth referred to as muon 
candidates), and re quir e that at least one is positively 
identified as a muon [16| . We require that the muon can- 
didates form a geometric vertex with Xvtx < 20 for 1 
degree of freedom and displaced transversely by at most 
2 cm [l3| from the nominal location of the e+e" interac- 
tion region. We perform a kinematic fit to the T can- 
didate formed from the two muon candidates and the 
energetic photon. The GM energy of the T candidate is 
constrained, within the beam energy spread, to the to- 
tal beam energy and the decay vertex of the T is 
constrained to the beam interaction region. We select 
events with —0.2 < ^/s — m{T) < 0.6 GeV and place a 
requirement on the kinematic fit Xr ^ (^^r 6 degrees 
of freedom) . We further require that the momenta of the 
dimuon candidate A^ and the photon are back-to-back in 
the GM frame to within 0.07 rad, and that the cosine of 
the angle between the muon direction and A^ direction in 
the center of mass of the A^ is less than 0.92. The selec- 
tion criteria are chosen to maximize e/VS, where e is the 
average selection efficiency for a broad m^" range and B 
is the background yield in the off-resonance sample. 

The criteria above select 387,546 T{2S) and 724,551 
T(3S') events (mass spectra for T{2S) and T(3S') datasets 
are shown in Fig. [3] in 18[). The backgrounds are dom- 
inated by two types of QED processes: "continuum" 
e+e^ — !■ 7/i+/i^ and the initial-state radiation (ISR) pro- 
duction of p*^, (/), J /"ip, '0(25'), and T(15) vector mesons. 
In order to suppress contributions from the ISR-produced 
— > TT+TT^ final state in which a pion is misidentified 
as a muon (probability ^3%/pion), we require that both 
tracks are positively identified as muons when we search 
for A^ candidates in the range 0.5 < to^o < 1.05 GcV. 
Finally, when selecting candidate events in the rjb re- 
gion with dimuon invariant mass to^^ ~ 9.39 GeV in 
the T(25) (T(35)) dataset, we suppress the decay chain 
T(25) ^ 72X6(1^, XbilP) -> 71^^(15) (r(35) ^ 
72X6(2-P), X6(2P) 7iT(15)) by requiring that no sec- 
ondary photon 72 above a GM energy of E2 =0.1 GeV 
(0.08 GeV) is present in the event. 

We use signal Monte Garlo (MG) samples [H, ^ 
T{2S) jA° and T(35) ~* jA° generated at 20 val- 
ues of TO^o over a broad range 0.212 < m^o < 9.5 GeV 
to measure the selection efficiency for the signal events. 
The efficiency varies between 24-55%, depending m^o. 

We extract the yield of signal events as a function of 
TO^o in the interval 0.212 < m^o < 9.3 GeV by per- 
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forming a series of unbinned extended maximum likeli- 
hood fits to the distribution of the reduced mass mji = 
'fjm'^if^ — 4rnj^. The likelihood function contains contri- 
butions from signal, continuum background, and, where 
appropriate, peaking backgrounds, as described below. 
For 0.212 < m^o < 0.5 GeV, we fit over a fixed in- 
terval 0.01 < mji < 0.55 GeV; near the J/ip reso- 
nance, we fit over the interval 2.7 < mR < 3.5 GeV; 
and near the ^1){2S) resonance we fit over the range 
3.35 < mil < 4.1 GeV. Elsewhere, we use sliding inter- 
vals fj, — 0.2 < mji < ^ -I- 0.1 GeV, where fi is the mean of 
the signal distribution of mR. We search for A" in fine 
mass steps Am^o = 2-5 MeV. We sample a total of 1951 
771^0 values. For each m^jo value, we determine the BF 
products Bns = B{r{nS) 7^°) x S^^, where n = 2, 3. 
Both the fitting procedure and the event selection were 
developed and tested using MC and T(4S') samples prior 
to their application to the T{2S) and T{3S) data sets. 

The signal probability density function (PDF) is de- 
scribed by a sum of two Crystal Ball functions 22] with 
tail parameters on either side of the maximum. The sig- 
nal PDFs are centered around the expected values of 
mn and have a typical resolution of 2-10 MeV, which 
increases monotonically with m^o . We determine the 
PDF as a fmiction of m^o using the signal MC samples, 
and we interpolate PDF parameters and signal efficiency 
values linearly between the simulated points. We deter- 
mine the uncertainty in the PDF parameters by compar- 
ing the distributions of the simulated and reconstructed 
e+e" JisRJ/i^, J/'>P ^ events. 

We describe the continuum background below mj^ < 
0.23 GeV with a threshold function /bkg("7fl) oc 

tanh (^2i=i Pi^Fi^ ■ parameters are fixed to the 
values determined from the fits to the e+e^ jfi'^^^ 
MC sample [ilj and agree, within statistics, with those 
determined by fitting the r(25), T(35'), and T{AS) sam- 
ples with the signal contribution set to zero. Elsewhere 
the background is well described in each limited mji 
range by a first-order {mn < 9.3 GeV) or a second-order 
{nin > 9.3 GeV) polynomial with coefficients determined 
by the fit. 

Events due to known resonances (j), J/V', '0(25'), and 
T{1S) are present in our sample in specific mn intervals, 
and constitute peaking backgrounds. We include these 
contributions in the fit where appropriate, and describe 
the shape of the resonances using the same functional 
form as for the signal, a sum of two Crystal Ball func- 
tions, with parameters determined from fits to the com- 
bined T{2S) and T{3S) dataset. The contribution to the 
event yield from </> — > K~^K^ , in which one of the kaons 
is misidentified as a muon, is fixed to 111 ± 24 (T(25)) 
and 198 ± 42 {T{3S)). We determine this contribution 
from the event yield of e+e^ 7^, (p K^K^ in a 
sample where both kaons are positively identified, cor- 
rected for the measured misidentification rate of kaons 



as muons. We do not search for candidates in the 
immediate vicinity of J/ip and ip{2S), excluding regions 
of ±40 MeV around J/ip ±5a) and ±25 MeV (« ±3cr) 
around ip{2S). 

We compare the overall selection efficiency between 
the data and the MC simulation by measuring the ab- 
solute cross section da/drriR for the radiative QED pro- 
cess e+e~ — > 7^"'"/i~ over the broad kinematic range 
< mji < 9.6 GeV, using the off-resonance sample. We 
use the ratio of measured to expected cross sections 
to correct the signal selection efficiency as a function of 
TO^o. This correction ranges between 4-10%, with a sys- 
tematic uncertainty of 5%. This uncertainty accounts 
for effects of selection, reconstruction (for both charged 
tracks and the photon), and trigger efficiencies. 

We determine the uncertainty in the signal and peak- 
ing background PDFs by comparing the distributions of 
w 4000 data and MC e+e" 1iskJ/'4>, J/i^ ^ M^M" 
events. We correct for the observed difference in the 
width of the mn distribution (5.3 MeV in MC versus 
6.6 MeV in the data) and use half of the correction to 
estimate the systematic uncertainty on the signal yield. 
This is the dominant systematic uncertainty on the sig- 
nal yield for m^o > 0.4 GeV. We estimate that the 
uncertainties in the tail parameters of the Crystal Ball 
PDF contribute less than 1% to the uncertainty in signal 
yield based on the observed variations in the J/t/' yield. 
The systematic uncertainties due to the fixed continuum 
background PDF for < 0.23 and the fixed contri- 
bution from e+e^ 70 do not exceed (Tbkg(i5ns) = 
0.2 X 10~^. These are the largest systematic contribu- 
tions for 0.212 < ruAo < 0.4 GeV. 

We test for possible bias in the fitted value of the signal 
yield with a large ensemble of pseudo-experiments. The 
bias is consistent with zero for all values of m^o , and we 
assign a BF uncertainty of (Tbias(Sris) = 0.05 x 10^^ at 
all values of m^o to cover the statistical variations in the 
results of the test. 

To estimate the significance of any positive fluctua- 
tion, we compute the likelihood ratio variable 5(to^o) = 
sign(7Vsig) 2 log(Lmax/ -^0) , where L^ax is the maximum 
likelihood value for a fit with a free signal yield centered 
at TO^o, TVsig is that fitted signal yield, and Lq is the 
value of the likelihood for the signal yield fixed at zero. 
Under the null hypothesis S is expected to be normal- 
distributed with /X = and a = 1 (Fig. [T]). Including 
systematics, the largest S values are 3.1 {T{2S)) and 2.8 
{T{3S)), consistent with a null-hypothesis distribution 
for 1951 TO^o points. 

Since we do not observe a significant excess of events 
above the background in the range 0.212 < m^o < 
9.3 GeV, we set upper limits on B2S and 633- We add sta- 
tistical and systematic uncertainties in quadrature. The 
90% confidence level (C.L.) Bayesian upper limits, com- 
puted with a uniform prior and assuming a Gaussian like- 
lihood function, are shown in Fig. [5] as a function of mass 
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FIG. 1: Distribution of the log-likelihood variable 5 with both 
statistical and systematic uncertainties included for (a) T{2S) 
fit, (b) r{3S) fit, and (c) combination of r{2S) and r(35') 
data. There are no points outside of displayed region of S. 
The solid curve is the standard normal distribution. 



ruAo . The limits vary from 0.26 x 10"^ to 8.3 x 10"^ (62s) 
and from 0.27 x lO^^ to 5.5 x 10"^ (63s)- 

The BFs B{r{nS) -fA°) are related to the effective 
coupling fr of the bound 6-quark to the through 0, 



B{T{nS) 7^") 
B{r{nS) 1+1-) 




'r(nS') 



(1) 



where I = e ov ji and a is a fine structure constant. 
The effective coupling fr includes the Yukawa coupling 
of the &-quark and the m^o -dependent QCD and rela- 
tivistic corrections to BnS [ll] and the leptonic width 
of T(nS) To first order in as, the corrections 

range from to 30% but have comparable uncer- 
tainties [2^. The ratio of corrections for T{2S) and 
T{iS) is within 4% of unity [2^ in the relevant range 
of m^o. We do not attempt to factorize these contribu- 
tions, but instead compute the experimentally- accessible 
quantity f^B^^ and average T{2S) and T{'iS) results, 
taking into account both correlated and uncorrelated un- 
certainties. The combined upper limits are shown as a 
function of m^o in Fig. [2fc) (plots with expanded mass 
scales in three ranges of rn^o are available in Fig. |4][6] 
in [13]) and span the range (0.44 - 44) x 10"^, at 90% 
C.L. The combined likelihood variable {S) = {w2sS2S + 
wssSss) / \/ + w'^s shown in Fig. [TJ;, where WnS is 
the statistical weight of the T{nS) dataset in the aver- 
age. The largest fluctuation is (S) = 3.3. Our set of 1951 
overlapping fit regions corresponds to w 1500 indepen- 
dent measurements [2^. We determine the probability 
to observe a fluctuation of (S) = 3.3 or larger in such a 
sample to be at least 45%. 

We do not observe any significant signal at m^o — 
0.214 GeV (Fig. [7] in [l3|) and set an upper limit on the 
coupling f^niAo = 0.214 GeV) < 1.6 x 10^^ at 90% C.L 
(assuming = 1), which is significantly smaller than 
the value required to explain the HyperCP events as light 
Higgs production [llj. 

A fit to the rji, region (Fig. [5] in |18j ) includes back- 
ground contributions from the ISR process e"*'e~ — > 




(GeV) 



B^^, (b) e(r(35) 



FIG. 2: 90% C.L. upper limits on (a) B{r{2S) jA°) x 
7^*^) X B^f^, and (c) effective coupling 
nij, aa a luiii^uiuii ui /it^u. Thc shaded areas show the 
regions around the J/^ and ^l'{2S) resonances excluded from 
the search. 



fx X Bfi^ as a function of m^o. 



liSR^i^S), and from the cascade decays T{nS) 
l2XbJ, XbJ ~> 7ir(lS') with T{\S) ■ We mea- 

sure the rate of the ISR events in the T(4S') dataset, 
scale it to the T{2S) and T(3S') data, and fix this con- 
tribution in the fit. The rate of the cascade decays, 
the number of signal events, and the continuum back- 
ground are free in the fits to the T{2S) and T(35) data 
sets. We measure B{T{2S) ~> 777^) x ^(77^ M^A^ ) = 
(-0.4 ± 3.9 ± 1.4) X 10-6 and S(T(3S') ^ ^Ib) x B{T]b 
fi+^j.-) = (-1.5 ± 2.9 ± 1.6) X 10-^ where the first 
uncertainty is statistical and thc second is systematic, 
dominated by the uncertainty in r(?7f,). Taking into 
account the BaBAR measurements of B{T{2S) ^-qb) 
and B{T{iS) ^ 777b) [3], we derive ^(776 ^ M+M") = 
(-0.25 ± 0.51 ± 0.33)% and ^(77^ ^ /x+zx") < 0.9% at 
90% C.L. This limit is consistent with the mesonic inter- 
pretation of the 77b state. 

In summary, we find no evidence for the dimuon decays 
of a light scalar particle in radiative decays of T{2S) and 
T(3S') mesons. We set upper limits on the coupling f^ x 
Bf,f, for 0.212 < 771^0 < 9.3 GeV. Assuming « 1 
in the mass range 2777^ < 777^0 < 1 GeV, our results 
limit the coupling fr to be at most 12% of the Standard 
Model coupling of the b quark to the Higgs boson. Our 
limits rule out much of the parameter space allowed by 
the light Higgs [9] and axion [6] models. We also set an 
upper limit on the dimuon branching fraction of the 77;,. 
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APPENDIX: EPAPS MATERIAL 

The following includes supplementary material for the Electronic Physics Auxiliary Publication Service. 
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FIG. 3: Distribution of the reduced mass m_R in (a) the T{2S) data and (b) the T{3S) data. Blue (open) liistogram shows 
the distribution for the selection in wliich only one of two muons is required to be positively identified. The peak from 
e'^e~ — > 7iSH,p''(770), p" — > Tv'^n~, in which one of the pions is misidentified as a muon, is clearly visible. Black (filled) 
histogram shows the distribution for the selection in which both muons are positively identified (this selection is used in the 
search for 0.5 < m^o < 1.05 GeV). The ISR-produced peaks J/-tp fJ-^ f^" and T(1S) — > p^/i" are visible for both selections. 
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FIG. 4: 90% C.L. upper limits on the effective Yukawa coupling fx x as a function of m^o in the range 0.212 < m^o < 
1.05 GeV. 
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FIG. 5: Upper limits on the effective Yukawa coupling fy x as a function of m^o in the range 1 < m^o < 4 GeV. The 
shaded areas show the regions around the J/tp and tp{2S) resonances excluded from the search. 
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FIG. 6: Upper limits on the effective Yukawa coupling fy- x Bfj,^ as a function of m^o in the range 4 < rrij^o < 9.3 GeV. 
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FIG. 7: Fits to the HyperCP region m^o = 0.214 GeV in (a) T{2S) dataset and (b) T{3S) dataset. The bottom graph shows 
the mu distribution (solid points), overlaid by the full fit (solid blue line). Also shown are the contributions from the signal at 
m^o = 0.214 GeV (solid red line) and the continuum background (dashed black line). The inset zooms in on the signal region. 
The top plot shows the normalized residuals p — (data — fit)/o-(data) with unit error bars. The individual fits correspond to 
the log-likelihood ratios of iS = 2.0 {T{2S)) and S — 0.2 (T{3S)), and the combined significance (<S) — 1.4. 
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FIG. 8: Fits to the rib region in (a) T{2S) dataset and (b) T(3S) dataset. The bottom graph shows the rriR distribution (solid 
points), overlaid by the full fit (solid blue line). Also shown are the contributions from the signal at at rrirj^ = 9.389 GeV (solid red 
hue), background from the e'^e" 7isrT(1S) (dot-dashed green line), background from T(3S) — > 7X6(2P), Xb{^P) — > I'^i^S) 
(dotted magenta line), and the continuum background (dashed black line). The inset zooms in on the signal region. The top 
plot shows the normalized residuals p = (data — fit)/a(data) with unit error bars. 



